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Abstract 

We investigate the cosmology of a recently proposed deformation of Einstein gravity, emerging 
from quantum gravity heuristics. The theory is constructed to have de Sitter space as a vacuum 
solution, and thus to be relevant to the accelerating universe. However, this solution turns out 
to be unstable, and the true phase space of cosmological solutions is significantly more complex, 
displaying two late-time power-law attractors - one accelerating and the other dramatically de- 
celerating. It is also shown that non-accelerating cosmologies sit on a separatrix between the two 
basins of attraction of these attractors. Hence it is impossible to pass from a decelerating cosmol- 
ogy to an accelerating one, as required in standard cosmology for consistency with nucleosynthesis 
and structure formation and compatibility with the data inferred from supernovae la. We point 
out that alternative models of the early universe, such as the one investigated here might provide 
possible ways to circumvent these requirements. 
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I. INTRODUCTION 



The observed acceleration of the universe, confirmed from a cross-correlation of indepen- 
dent cosmological datasets a crucial question for fundamental physics. 
There are at least three distinct approaches from which an explanation of these observations 
might come. The simplest explanation is that cosmic acceleration is due to a cosmological 
constant. If this is the case, then our understanding of quantum fluctuations of matter fields 
and their gravitational effects will require new insights. A second possibility is that there is 
no cosmological term and that the dynamics of some new matter field, for example a scalar, 
lead to late-time acceleration. Such an explanation would require us to understand the exis- 
tence of an unnaturally weakly coupled long-range field, with minute energy density at this 
particular epoch in cosmic history. A third, and less well-explored, possibility is that cosmic 
acceleration is due to an infrared modification of gravity, which yields cosmological solu- 
tions radically different from those of pu re g eneral relativity. While many such models exist 

it has recently been shown by two of the authors that a unique theory is selected by the 



special geometry of Lorentzian spacetimes approximating a quantum spacetime |33l . 13^ 

In this paper, we explore the cosmology of this distinguished deformation of Einstein 
gravity. The deformation takes the form of an inverse curvature term constructed from 
the Riemann tensor, similar to, but distinct from the modifications discussed in [9|. In 
particular we investigate in some detail the phase space structure of the resulting modified 
Friedmann equation. In ordinary Einstein gravity, if the present acceleration of the Universe 
is attributed to a cosmological constant, then any cosmology inevitably approaches the de 
Sitter vacuum at late times. In our case, although the theory is constructed to have de 
Sitter space as a vacuum solution, and thus to be relevant to the accelerating universe, this 
solution is unstable, and the true phase space of cosmological solutions is significantly more 
complex, displaying two late-time power-law attractors - one accelerating and the other dra- 
matically decelerating. However, it appears that the theory cannot explain the accelerating 
universe. This is because non-accelerating cosmologies sit on a separatrix between the basins 
of attraction of the two attractors. It is therefore impossible to pass from the decelerating 
cosmology required, in standard cosmology, by nucleosynthesis and structure formation to 
the accelerating one inferred from supernovae la. The supernovae data provides evidence 



2 



for the jerk at greater than 99% confidence level 35]. The phase space is divided into ever 
accelerating and ever decelerating universes; cross-overs cannot occur. However, for alter- 
native cosmologies, such as the one under investigation here, the crossover requirement is 
less clear and there might be some loopholes which we will discuss. 

The structure of this paper is as follows. In the next section we describe inverse Riemann 
gravity and then derive the cosmological equations of motion necessary for the remainder 
of the paper. In section IIHI we provide an analytic phase space description of the most 
important features of the solutions. We show that the de Sitter solution is unstable, identify 
late-time power-law attractors and, most importantly, demonstrate the existence of a critical 
zero-acceleration separatrix. In section HVI we go on to provide complete numerical solutions 
to the equations of motion, mapping out the phase space for both the vacuum and non- 
vacuum cases, before concluding with a discussion in section 



II. INVERSE RIEMANN COSMOLOGY 



A. Inverse Riemann gravity 

Quantum gravity heuristics suggest that Lorentzian manifolds approximating a quantum 
spacetime possess both infrared and ultraviolet sectional curvature bounds [3J2, yj]. Such 
bounds are entirely equivalent to bounds on the eigenvalues of the Riemann tensor, consid- 
ered as an endomorphism on the space of 2-forms. If one attempts to write down an action 
such that all solutions respect these bounds, one quickly discovers that there is a large class 
of such theories (as many as there are holomorphic functions on an annulus). The ambiguity 
arises because of our ignorance of the exact quantum spacetime structure in this approach. 

However, because we are interested in cosmology, it is appropriate to take the classical 
limit, corresponding to the removal of the curvature bounds. If one further requires that, 
in this limit, de Sitter spacetime is a solution, as one might hope occurs in order to explain 
the accelerating universe, then the resulting classical theory is unique j^j]. Surprisingly, 
and interestingly, it is not the Einstein-Hilbert action with a positive cosmological constant. 
In this paper, we take the resulting classical action as our starting point, and proceed to 
investigate its broad cosmological characteristics. 

The action is a new type of infrared modification of Einstein-Hilbert gravity. Including 
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a matter action S m , the total action in d spacetime dimensions reads 

S = jd d x^TT^R- ^C'iT 1 ) + S m , (1) 

where R b is the endomorphism on the space of antisymmetric two-tensors defined by 
the Riemann tensor R^- ab \ c ^. The trace should not double-count, and hence, is defined by 
Trf(R) = f(R) A A = /(-R)' a ^[ ab ]/2. We have chosen units such that the reduced Planck mass 
is Mp = (87rG) _1//2 = 1. The vacuum solutions are de Sitter and anti de Sitter space of 
constant curvature, proportional to ±( respectively. 

The action ((TJ is similar to that introduced in [?J]. However, note that the current 
action contains both the Ricci scalar and a specific function of the full Riemann tensor, 
rather than just a function of the Ricci scalar. In particular, it is important to note that 
Tr(i? _1 ) is generally not equal to the inverse of the Ricci scalar. As noted in j^, the 
dimensionality of the space of metrics is different from that of the space of Riemann tensors, 
which means that it is not possible to perform a transformation of the metric, such as a 
conformal transformation, which reduces the action to that of Einstein gravity, with the 
extra degrees of freedom represented by additional matter fields. This means that direct 
comparison with solar system tests is a more complicated calculational task, which we will 
not attempt in this paper. It thus remains to be seen whether this action is compatible with 
gravitational dynamics at the solar system and galactic scales. 

The general equation of motion resulting from (0) is 



1 d — 2 

T d + 2 s 



(iri)<w) 6 + x -g^{R^y\ h - v fe v c (/r 2 ) c fe> 



T ij , (2) 



where T u is the matter energy-momentum tensor derived from S m , and the sign convention 
R a bcd = d c T a bd + r a ec r e M - (c d) has been used. 

B. Cosmological equations of motion 

We impose the cosmological Friedmann-Lemaitre-Robertson- Walker ansatz 

ds 2 = -dt 2 + a(t) 2 rfS 2 (3) 

for a homogenous and isotropic spacetime, where a(t) is the scale factor and 

dv 2 

dZ 2 = g a pdx a dx p = — + r 2 (d9 2 + sin 2 9 dcj) 2 ) (4) 

1 rhV 



is the metric on a d — 1 dimensional space of normalized constant curvature k = 0, ±1 with 
Riemann tensor Rap^s = k(g ai gps — g a sg^)- With this ansatz, R a b is a diagonal matrix, 
with indices taking values in {[0/3], [a/3]}. The corresponding d(d — l)/2 eigenvalues are 
given by 



R m m 

R [a % s] 



V 

a" 
(k + a 2 



(5) 
(6) 



the overdot denoting differentiation with respect to cosmic time t. Thus, the inverse curva- 
ture correction to the Einstein-Hilbert action appearing in is given by 



Tr (iT 1 ) = (d - 1) 



d-2 



k + as 



a 

+ - 



(7) 



From the form of this correction we may already anticipate a possible difficulty in passing 
from a decelerating phase to an accelerating one, since (JJJ) typically blows up when a — 0. 
We consider perfect fluid matter with the energy-momentum tensor 



T ij = (p + p)«V + pg l 



(8) 



where u % characterizes the fluid's four-velocity in its rest frame, p denotes its energy density, 
and p its pressure. The analogue of the Friedmann equation is derived from the time-time 
component of the equations of motion which, after some algebra, becomes 



2C 2 
d+2 
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(9) 



a 2 (d-l)(d-2) 
Note that this result, from the full equations of motion, agrees with an effective action 
calculation performed with the line element (J3J) additionally including a lapse function N(t) 
which is set to unity after variation. Setting £ = 0, one obtains the usual Friedmann equation 
of the standard cosmology based on the Einstein-Hilbert action. For later convenience, we 
adopt the usual definition of the Hubble parameter H(t) = a(t)/a(t), which allows us to 
rewrite the above equation in the form 



2C_ 
d + 2 



H + H 2 



+ d 



H 



H + H 2 
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(10) 



This equation and the continuity equation 

p+(d-l)H(p + p) = 0, (11) 

which follows directly from the conservation of energy-momentum V a T a6 = 0, are equivalent 
to the full set of gravitational field equations. 

III. PHASE SPACE DESCRIPTION OF SOLUTIONS 

The main goal of this work is to understand the expansion history of the universe in 
the theory introduced in the previous section. Because the resulting system of equations 
is of order greater than two, it is most convenient to understand the generic properties by 
considering its phase space evolution. 

As is typical in cosmology, we close the system of equations by specifying a matter 
equation of state p = p(p). Then, without further assumptions, the cosmological phase 
space of the resulting system is given by (p, a, a, d). The diffeomorphism invariance of the 
field equations (in particular, the time reparametrization invariance) bundles the trajectories 
in this phase space into classes of physically indistinguishable solutions. We may pick one 
member of each class by requiring that for a fixed time to the scale factor should satisfy 
a(t ) = 1. Note that this does not change the dimension of the phase space; it merely 
implies that, in finding solutions, one does not need to scan over different initial conditions 
for the scale factor a. 

Note that upon a global conformal rescaling \— > £ _1 <7y, the Riemann endomorphism 
scales with £ and the metric density with £~ d / 2 . This allows to factor out £( 2_d )/ 2 while 
simultaneously rescaling the deformation parameter ( — > For a matter action that also 
scales with £( 2_d )/ 2 ; such as a massless scalar or simply the vacuum, solutions of theories 
with different non-zero ( are therefore related by global conformal transformations. This 
means that the qualitative nature of the phase space will be the same for all values of £. 

We now survey the salient features of the phase space. 

A. De Sitter vacuum instability 

The equations of motion (J2J) have precisely two solutions of constant curvature, which 
are the de Sitter and anti-de Sitter spacetimes with Hubble parameter (in the flat slicing) 
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±£. Note that this statement is independent of the spacetime dimension d. 

For convenience we define J = H and solve equation (fTU|) with k = for H. This yields 
the system of equations 



P 
H 

j 



-(d-l)Hp, 
J, 

(J + H 2 f \ d + 2 
2C 2 



+ 



2H 
dH 2 



H 2 



•2f> 



(d- l)(d-2) 



J + H 2 



2H(3HJ + H 3 ) 3(d-2) 



(12) 
(13) 

(14) 



(J + H 2 ) 2 (J + H 2 ) 3 2H 2 

This system possesses two fixed points. The first is at (p, H, J) = (0, \/C, 0), and coincides 
precisely with the de Sitter vacuum solution a(£) = exp y/(t. The second fixed point is 
(p, iJ, J) = (po, 0, 0) for constant p . However, this fixed point satisfies J + H 2 = 0, which 
means that it occurs at a singularity of equation ([10)1 . Therefore this is not an admissible 
solution. 

To analyze the stability properties of the de Sitter fixed point under cosmological evo- 
lution, we linearize the above equations setting (p, H, J) = (x, v^C + V-> z ) with small 
X T = (x,y,z). As we shall see, the de Sitter solution is a saddle point in phase space, 
and hence unstable. 

The linearized system can be written in the form X = MX with the matrix M given by 

(-(d-i)VC o o \ 

1 

(d + 2)( -(d-i)vT 



M 



\ 



(15) 



/ 



2{d-\)(d-2) 

The eigenvalues of this matrix determine the type of fixed point; they are 



Ai 

A± 



\Vl(-d+ 1± \f¥+2d + V 



(16) 
(17) 



Clearly, Ai < 0, implying the stability of the de Sitter solution against perturbations of the 
energy density p away from the vacuum value. However, A + > and A_ < 0, implying that 
this solution is a saddle point, having one unstable direction in phase space. This behaviour 
can be seen clearly in the numerical phase space plots presented below. 
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B. Power law attractors 



The simplest form of accelerating expansion in the final stage of the cosmological evolu- 
tion, namely a stable de Sitter vacuum, is ruled out, as we have seen. However, there exist 
power law attractors. 

Substituting the ansatz a ~ t 7 , with constant 7, into the modified Friedmann equation 
(JHJ), one finds late-time solutions (t — > 00). For k = 0, in d > 2, we have 

7 ± = ^(3rf±v / 6(rf 2 + 2)). (18) 

In d = 4, these exponents are 

7 i d=4) =2± v / 3, (19) 
and correspond to effective equation of state parameters 



(20) 



cS ~ 3 " 

We see that 7+ > 1, which describes an accelerating solution, and < 7^ < 1, which 
describes a decelerating solution. For k 7^ 0, there is only the accelerating late-time solution 
7 + , since the curvature term ultimately dominates, although, if inflation vastly flattens 
space at early times, then the decelerating behavior during the rest of cosmic history may 
include an intermediate regime approximated by the decelerating solution before curvature 
domination in the far future. A simple linearized analysis about these solutions shows that 
they are both attractors, with perturbations from them decaying away at late times. 

The saddle point instability of de Sitter is directly connected to the existence of the 
accelerating power law attractor. As we will see, a large class of cosmological phase space 
trajectories will come close to the de Sitter vacuum for a time, but ultimately approach 
power law behavior. 



C. The acceleration separatrix 

Perhaps the most cosmologically relevant feature of the modified Friedmann equation Q 
is the generic singularity at a = 0. This is, in fact, the location of a phase space separatrix. 
To see this, note that near a = 0, equation (JHJ) takes the asymptotic form 

j t a{t) ~ 2H(t)d(t) . (21) 
8 



Since the Hubble parameter is always positive, it is clear that already accelerating solutions 
(with a small and positive) are driven to accelerate even more, and already decelerating ones 
(with a small and negative) are driven to further deceleration. Thus, phase-space trajectories 
can never cross over from deceleration to acceleration. 

This seems to create an almost insurmountable obstacle to the construction of a realistic 
cosmology in which the late-time acceleration of the universe is explained by the modified 
gravity nature of this theory. In standard cosmology, successful nucleosynthesis and struc- 
ture formation require a decelerating cosmology at early times and hence a transition from 
deceleration to acceleration is required by observations. 

It is worth mentioning, however, that much of the intuition we have from standard cos- 
mology is based on the implicit assumption that the Einstein equations (via the unmodified 
Friedmann equation) control the background dynamics. If instead, a modified Friedmann 
equation is used, many calculations beyond the scope of this paper must be redone in order 
to conclusively rule out a modified gravity theory. These include structure formation simu- 
lations and the determination of cosmological parameters, such as the total matter density, 
from a fit of the modified model to the cosmic microwave background. One clear example in 
which the standard wisdom proves inapplicable is the following: successful nucleosynthesis 
is possible even in an accelerating background if there is an asymmetry in the neutrino and 
antineutrino abundances. This is shown in 35], without assuming a specific background 
dynamics. 



IV. NUMERICAL SOLUTIONS 



While the asymptotic behavior of the theory can easily be analyzed analytically, we must 
solve the full equations numerically. To get a feeling for how solutions look, we restrict our 
attention to the case of d = 4, flat k = cosmologies which have the most interesting phase 
space features. We first plot solutions to the vacuum equations of motion, setting p = p = 0. 
Figure H] shows solutions in the (H, H) phase space. In the phase space plot we have scaled 
( so that the unstable de Sitter solution is located at the point (0,1). The two late-time 
power law attractors and the separatrix corresponding to a = are clearly visible in the 
upper left quadrant of the phase space. Solutions with power-law behavior (a ~ t 7 ), follow 
the curves H = y —^H in the phase space. The upper right quadrant of the phase space is 
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included for completeness, however, these solutions have H > and are physically irrelevant 
for the purposes of our paper. 




dH/dt 



FIG. 1: Vacuum solutions for £ = Hq = 1. The de Sitter vacuum sits on the H = axis at this 
value (large dot). The two dashed curves mark the power law attractors. The accelerator is the 
upper of the two curves and the decelerator is the lower one. The separatrix is indicated by the 
dotted line. 

Since we are interested in the late universe, we now specialize to the case in which the 
dominant matter component is pure dust, for which the pressure vanishes p = 0. We set 
= 0.28, a = 1 and H(t = t ) = H , where a subscript denotes the value of a quantity 
today. Figure 121 shows these solutions in the same (H,H) phase space. 

We focus on solutions with initial conditions p(t Q ) and a(t Q ), obtained from observations 
of the current values of the Hubble parameter and the density parameter Q m = H~ 2 p/ (d—1), 
consistent with the present state of the universe within the context of Einstein gravity. We 
therefore set p(t ) = ?>HlVL m) a(t Q ) = 1 and d(t ) = H . It is then sufficient to scan over 
initial conditions for the remaining phase space variable a. Since we are in four dimensions, 
we only exhibit projections to three-dimensional hypersurfaces of constant a = 1. For the 
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H 




FIG. 2: Solutions to equation @, with Q$ = 0.28, for a variety of different initial conditions. 
The thick solid line denotes initial conditions that are consistent with constraints on the evolution 
of the scale factor from observations of type la supernovae. These constraints exclude the other 
solutions at the 95% confidence level. 



present values of the matter density and the Hubble parameter we choose Q m = 0.28 and 
H — 1. To find initial conditions for a (or equivalently H) compatible with the recent 
supernovae data, we assume Einstein gravity and a flat universe filled with dust (fl m ) and 
a "dark energy" component with effective instantaneous equation of state parameter u> - 
From the Einstein equations we find 



H = -H 1 + , (22) 

2 V 1+a ' 



where 

n = 

1-Q 



a = -^-. (23) 



7/1 



The initial conditions for H, compatible with supernovae data are located in the green line 
in Figure 2, and are obtained using (1221) together with the observational data —1.48 < < 



-0.72 at the 95% confidence level 



01 



37 
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V. CONCLUSIONS 



The discovery of cosmic acceleration seems to imply that there is something important 
that we do not understand about particle physics, gravity or the interaction between them. 
While there have been many attempts to understand this phenomenon at a phenomenological 
level, it would be extremely exciting to derive the late-time evolution of the universe from 
a fundamental theory. 

In this paper we have investigated the cosmological evolution of the classical modified 
gravity theory given by the action (0). This theory arises in the following way: considerations 
of quantum gravity heuristics suggest that Lorentzian manifolds approximating a quantum 



33|, |3J. A large class of dynamics whose 



spacetime possess sectional curvature bounds 
solutions respect these curvature bounds exists. For a certain subclass, the classical limit 
corresponding to the removal of the curvature bounds, leads to an interesting result: one 
obtains a unique theory containing de Sitter space as a vacuum solution. The action for this 
theory is not the Einstein-Hilbert action with a positive cosmological constant, but rather 
contains a correction in terms of the inverse eigenvalues of the Riemann tensor. 

Using the resulting action as our starting point, we have investigated its broad cosmo- 
logical characteristics. As is typical of modified gravity theories with inverse curvature 
terms 0, , the equations of motion contain terms that are singular for certain values of 
the metric components. In the cosmological setting this translates into certain combina- 
tions of time derivatives of the scale factor. In many models [j], Q] these singularities lie 
in regions of the parameter space that are never reached by trajectories corresponding to 
realistic cosmological evolution. However, in the specific theory defined by the action ([[J, 
the relevant singularity of the equations of motion occurs precisely when d(t)=0. 

If one interprets the available cosmological data sets under the assumption that Einstein 
gravity controls the dynamics of the background, there is good evidence that the universe 
underwent a decelerating phase (needed for successful nucleosynthesis and structure forma- 
tion), and now a late-time accelerating phase (to fit, among other data, the observations of 
type la supernovae). Thus standard cosmology requires that the universe pass from deceler- 
ation to acceleration at some point in cosmic history. This is not allowed by the singularity 
structure of the theory investigated here. 

Assuming that, even with the modified Friedmann equation (JHI), nucleosynthesis and 
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structure formation can only take place during a decelerating phase, as in the standard Big- 
Bang model, the exclusion of cross-overs from decelerated to accelerated expansion presents 
an insurmountable obstacle to a viable cosmology. It is worth commenting, however, that in 
order to reliably confront alternative cosmologies with observations, one must either use data 
without assuming specific background dynamics, as is done e.g. in (35, 39 1, or redo standard 
calculations (such as extracting the cosmological parameters from data on the basis of the 
modified dynamics). This of course presents a formidable task and is beyond the scope of 
this paper. 
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